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Abstract 

In this paper, we give a positive answer to a problem posed by Nakagawa, Sakamoto and Yamamoto concerning a 
nonlinear equation with a fractional derivative. 
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. 1. Introduction 

In their overview paper concerning the mathematical analysis of fractional equations, Nakagawa, Sakamoto and 
Oh I Yamamoto [ 1 ] posed the problem concerning global solutions and blowing-up in a finite time of solutions to the 
equation 

^ ; ( c B%u(t) = -u(t){\ - w(f)), f>0, 

\ u(0) = mo, 

C3 ■ 

where L Dy is the Caputo derivative defined for g e C [0, T] by 
forO < a < 1. 

in 

^| ■ Let us recall, in the case a = 1, the results concerning solutions of ([1) : 

j | . - For < u(Q) < 1, the solution exists globally. Moreover, 

o: i 
(N . i«( f )i < -771 — ^ °> as f - * +o °- 

e f (l - M ) 



- For m(0) > 1, the solution can not exist globally. 

Here, we show that the same conclusions are valid for equation ([TJ. Moreover we analyse : 

1. The large time behavior of the global solution. 

2. The blow-up time and profile of the blowing-up solutions. 

Note that if we set w — u - 1, then <[TJ reads 

c D%w(t) = w(f)(l + w(tj), 

which describes the evolution of a certain species; the reaction term w(l + w) describes the law of increase of the 
species. 
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2. Preliminaries 

In this section, we present some definitions and results concerning fractional calculus that will be used in the 
sequel. For more information see 

The Riemann-Liouville fractional integral of order < a < 1 of the integrable function / : M + — > R is 

J if(ty-=^ [(t-rr-'mdr, t >o, 

Ha) Jo 

where T(a) is the Euler Gamma function. 

The Riemann-Liouville fractional derivative of an absolutely continous function fit) of order < a < 1 is 

D o + /« : = ^77—, i fit ~ rVm dr. 
T(l - a) dt Jo 

The Caputo fractional derivative of an absolutely continous function /(f) of order < a < 1 is defined by 

c D a 0+ f(t) := J^'f'it) = f7T~ 7 fit ~ T)- a f'iT) dr. 
Til -a) Jo 

Both derivatives present a drawback : 

- The Riemann-Liouville derivative of a constant is different from zero, 

D a C * 0, 

while the Caputo derivative require fit) to calculate C D" ) fit), for < a < 1. 

- We know that the Riemann-Liouville derivative of the Weierstrass function exists for any < a < 1, but not 
for a = I, 

But for regular function with /(0) = 0, both definitions coincide. 

Next, we recall a lemma that will be used hereafter. 

Lemma 2.1. ( see i3/J. Let a, b, K, i]/ be non-negative continuous functions on the interval I = (0, T) (0 < T < <x>), 
let a) : (0, oo) — > R be a continuous, non-negative and non-decreasing function with a»(0) = and a>iu) > for 
u > 0, and let Ait) — maxo< s <r a(s) and Bit) = maxo< s <r bis). Assume that 



•f- 

Jo 



4>it) < ait) + bit) K(s)a)(ifr(s)) ds, t e /. 
Jo 



Then 



HiAit)) + Bif) Kis)ds 



f. 

Jo 



, fe(0,7U 



where Hiv) = JT" ^ (v > v > 0), H 1 is the inverse ofH and T\ > is such that HiAit)) + Bit) £ Kis)ds e DiH 1 ) 
for all t € (0, tI). 

Here, we consider the problem 

~D a uit) = -u(t)(l - uit)), 



K(0) = Ho, (2) 



for < a < 1 and uq > 0. 
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3. Main results 

The local existence of solutions to © is assured by the 
Theorem 3.1. (see We consider the fractional differential equation of Caputo's type given by 

c D a Q+ u(t) = f(t,u(t)), t>0, 
u(0) = uq. 

For < a < I, u e R, b > and T > 0. 
Assume that 

1. / € C(R , M) where R = {(t, u), < t < T, \u - u Q \ < b] and \f(t, u)\ < M on R ; 

2. \f{t, u) - f{t, v)| < L\u - v|, L > 0, (?, u) € R . 

Then there exists a unique solution u e C([0, h])for (0, w/iere /z = min jr, ^ — — - j J. 
Theorem 3.2. Lef m foe the solution of problem d2j. We have : 

- If < uq < 1, f/ie solution is global and it satisfies < W < 1. Moreover, u is given by 

u(t) = E a (-t a )u + f (t- s) a - l E a , a (-(t - s) a ) u z (s)ds, 
Jo 

and for some constants c > and c\ > 0, we have 

<u(t)<- — , t> T :=(-?— )". 

-! ZLp \c\cuqI 

cuq a 

- If uq > 1, the solution blows-up in a finite time T* : lim u(t) — +oo. 
Moreover, we have the bilateral estimate : 

w(t) + 1 < u(t) < w(t) + I, 

and 



w/zere 

w(0 + i ~ ?r^(^ - ty a , as t^T~ x 

2 T(a) 

w(t) ~ ^^-(T- - ty a , as t — > T w . 
Here, T„ is the blow-up time of w, which satisfies 

r (a+i M°<r,<f^)°. 

U(« -i)i - H -\ u -\ ) 
and 7V is the blow-up time of w, which satisfies 

(I^-f <T,<( Y -^-t 
U(mo-1)/ " \ «o-l / 



(3) 
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Proof of Theorem 3.2. 

Part 1. If < uq < 1, then the solution is global. 
The solution to (f2j) is given by 



u(t) = E a (-t a )u Q + f (f-i) ff -'£ ff>ff (-(f-i) a )M z (i)£/i. (4) 
Jo 

Where the Mittag-Leffler functions E a (-f) and E a>a {— f) are defined by : 

™ {-\)ifi 

If m > 0, then u(i) > as £<*(-*") > and E a , a (-f) > 0. 

Now, we set the function 5(f) = 1, f > 0. 

As < mo < 1, then mo < 17(0). In addition, we have 

c Dlu(t) = = -5(0(1 - 5(f)). 
Hence 5 is an upper solution of the equation (0, and we have u(t) < 5(f) = 1, (see Q], Thm. 2.4.3, p. 32). 

Now, we examine the large time behavior of the global solution < u < 1 . 
For, let us recall the estimates ( see J5t]) : 

- For < a < 1, there exists a constant c > such that, 

< E a (-f) < — — <c, t > 0. (5) 
\+t° 

- For < a < 1, there exists a constant c\ > such that 

0<t a - x E a J-f)< Cl t a -\ f>0. (6) 
From © and using the inequalities (0 and ©, we obtain 



u(f) < cu +c\ f (f - 5) ff L u\s)ds. (7) 
Jo 



We apply Lemma 2.1 to © with co(x) = x 2 , K(s) = (f - s)" -1 , A(t) = cu , B(t) = ci. 
For t > To, we have 

H(cu ) + — f € D(H' 1 ), 

a 

where //(v) = and // _1 (z) = — , z # — . 

v v -±- - z v 

"0 

So we obtain, 



w(0 < # _1 



ff(CKo) + — f" 



Therefore 

M (f) < — , f > r . 

_! £i fa 

cuq a 

Part 2. If wo > 1, then the solution blows-up in a finite time. 



1 . We show that u > 1 . For, let us define the new unknown function w = u — 1 . The function w satisfies 

c D»w(t) = w(f)(l + wit)), 

w(0) : = Wo = Uo - 1. 

As wo > 1, then wo > 0. Moreover, we have (Q) 

w(?) = £ ff (f> + f (f - s) a - l E aia ((t - s) a ) w 2 (s)ds. 
Jo 

Therefore, w > 0; hence k > 1. 

2. We prove that u blows-up in a finite time. 

Since we have w(f) = «(f) - 1, it is seen that if u(f) — > oo as f — > 7"*, then w(f) — > oo as f — > T* and 
vice versa. That is w and m will have the same blow-up time. 

We now must examine the blow-up properties of w, the solution of problem ©. These are obtained by compar- 
ing wit) with the solutions of the following problems : 

c D«w(f) = w 2 (t), 



and 



w(0) = wo, 



c Dlw(t) = (w(t) + l)\ 
w(0) = wq. 



We see by comparaison (|01) that 



w(f) < wit) < w(t), < t < min{r F , 7~}. 

Following the paper of Kirk, Olmstead and Roberts 0, we may assert that the solution w ( resp. w) blows-up 
in a finite time TV ( resp. T„), such that 



\ 4wo / V wo / 



and 

So we have the following estimates 
Whereupon 



f r(a+1) )" < rs < (EffLtH) 4 

U(wo + i)/ Vwo + i/ 



T^<T* < Tyf. 



I Tia+Vrd < r < /IXa+l)\: 
V4(w + 5)/ ~ ~ \ w / 



□ 



4. Numerical implementation 

In this section, we will approximate the solution u given by ©. For, we need a numerical approximation of the 
convolution integral; this can be obtained using the convolution quadrature method. 
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Figure 1: solutions for a = 0.5 and «o = 5, 3, 2. 




Figure 2: solutions for uq = 5 and a = 0.3, 0.5. 



As it has been explained in [7], a convolution quadrature approximates the continuous convolution 

f Kit - s)f(s) ds, f>0, 
Jo 



by a discrete convolution with a step size h > 0. Then 



K(t„ - S)f(s)ds ~ ^ Un-jf(tj) 



where tj = jh, j = 0, 1, 2, n and the convolution quadrature weights u>j are determined from their generating power 
series as 

. £(£)j 

7=0 

Here : s} is the Laplace transform of K(t) and 5(0 is the generating polynomial for a linear multistep method. 

Let m„ be the approximation of u{t„) for n > 0. Using the convolution quadrature method we obtain 



M„ = (1 - W()) 1 £ , Q .(-f")M + 2^ W„_ ; M; 

/=0 



, « = 1,2,3.... 

Now, we introduce the following algorithm which gives the numerical approximation of solution to equation (0. 
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Algorithm 

Input : Give a, < a < 1 and uq, uq > 1. 

Initializations : Discretize the time with a step size h > 0; r,- = zTj, for all ;' = 1,2, n, z^pp* = Wo> «' = («o) 2 - 
Step 1 : Approximate the Mittag-Leffler function GML. 

Step 2 : Calculate convolution quadrature weights W using the fast Fourier transform (FFT). 

Step 3 : Calculate u' . 

do 

«* = GML *u l appx + Yi* «'■-'. 
«V = (1-W(1))- 1 *«*. 

«' = («Wx) 2 - 
i' = i' + l. 

until (m^x blows up) or (;' > n ). 

Output : Numerical approximation of w. 



Example 1. For Figurel, we set a — 0.5; f/ze initial conditions are respectively uq — 5, uq — 3 anc/ uq — 2. 

For Figurel, we take the initial condition uq — 5 and we plot the solutions; the dotted curve is the solution for a = 0.3 

and the solid curve corresponds to the solution for a = 0.5. 

As it has been proved, the solution blows up in a finite time which depends on Uq and a. 
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